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I. Introduction 


T he accurate calculation of aerodynamic forces and moments is of significant importance during the design 
phase of an aircraft. Reynolds-averaged Navier- Stokes (RANS) based Computational Fluid Dynamics 
(CFD) has been strongly developed over the last two decades regarding robustness, efficiency, and capabilities 
for aerodynamically complex configurations. 1,2 Incremental aerodynamic coefficients of different designs can 
be calculated with an acceptable reliability at the cruise design point of transonic aircraft for non-separated 
flows. But regarding absolute values as well as increments at off-design significant challenges still exist to 
compute aerodynamic data and the underlying flow physics with the accuracy required. 

In addition to drag, pitching moments are difficult to predict because small deviations of the pressure 
distributions, e.g. due to neglecting wing bending and twisting caused by the aerodynamic loads can result in 
large discrepancies compared to experimental data. Flow separations that start to develop at off-design con- 
ditions, e.g. in corner-flows, at trailing edges, or shock induced, can have a strong impact on the predictions 
of aerodynamic coefficients too. 

Based on these challenges faced by the CFD community a working group of the AIAA Applied Aerody- 
namics Technical Committee initiated in 2001 the CFD Drag Prediction Workshop (DPW) series resulting 
in five international workshops. The results of the participants and the committee are summarized in more 
than 120 papers. 3 7 The latest, fifth workshop took place in June 2012 in conjunction with the 30th AIAA 
Applied Aerodynamics Conference. 8 

All workshops were focused on the following key objectives: 6 

• assess state-of-the-art CFD methods as practical aerodynamic tools for the accurate prediction of 
forces and moments on industry-relevant aircraft configurations, with a focus on absolute as well as 
incremental values, 

• setup an international forum of experts from industry, research and academia for the verification and 
validation of RANS based CFD methods by applying different meshing methods and turbulence models, 

• define areas for additional research needed, 

• build, use, and maintain a public-domain transonic flow database for transport aircraft geometries 
including CAD data, grids, and numerical and experimental results, 

• document workshop findings and to disseminate through presentations and publications. 

NASA and the DLR Institute of Aerodynamics and Flow Technology are supporting these objectives as 
committee members and participants. 9-13 

The first three workshops used DLR transonic wind tunnel model configurations and experimental data 
achieved together with ONERA. 3 5,14,15 For the fourth and fifth workshops a new configuration the so- 
called Common Research Model (CRM) was defined by NASA and Boeing 16 for transonic flow conditions, 
see figure [fl In 2009/2010 experimental wind tunnel campaigns were performed with the CRM by NASA 
in the National Transonic Facility (NTF) at Langley and in the Ames 11-ft tunnel. The data have been 
published recently in several papers. 17-19 

A major aspect came into focus when the DPW-4 and DPW-5 computational results of the participants 
have been compared to the experimental data. Besides moderate discrepancies in drag at the cruise design 
point significant offsets of the pitching moments have been observed. These have been traced back to the 
model support system, which extends vertically from the aft fuselage, and to a deviation of wing twist 
between the computational and the wind tunnel model geometries. 20,21 

DLR results in DPW-4 and DPW-5 also showed differences between experimental data and numerical 
calculations of the wing pressure distributions, especially for the most outboard sections as presented in 
Figs. Hand ^nearly independent of the different grid refinement levels (L2: coarse, L6: ultra fine). Therefore, 
it is the first objective of these investigations to evaluate the influence of static aeroelastic wing deformations 
onto pressure distributions and overall aerodynamic coefficients. NASA and DLR decided to perform, in 
addition to their other investigations, fluid- structure-coupled simulations based on the NASA finite-element 
structural model of the CRM wind tunnel model and the DLR TAU CFD solver. 

A second aspect of DLR results in DPW-5 has been the prediction of a small flow separation (< 1% 
local chord length) at cruise conditions in the wing- fuselage junction near the wing trailing edge. The size of 
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Figure 1. NASA Common Research Model (CRM) without horizontal tail. 




(a) 77 = 0.501 


(b) 77 = 0.95 


Figure 2. Pressure distributions, TAU results for common grids (L2: coarse, L6: ultra fine), SA model, NASA 
NTF test data, Case 1. 




(a) 77 = 0.501 


(b) 77 = 0.95 


Figure 3. Pressure distributions, TAU results for common grids (L2: coarse, L6: ultra fine), Menter ku-SST 
model, NASA NTF test data, Case 1. 
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that separation and its development for increasing incidence angles was found to depend on grid resolution, 
topology, turbulence modelling, and numerical dissipation as demonstrated in Figs. H and Hi The effect is 
more pronounced for the Spalart-Allmaras turbulence model than for the Menter kcj-SST model. For very 
fine grid levels of both the hexahedral common grids and the DLR unstructured hex-dominant grids with 
overlapping hexahedral element corner block, separation starts close to the location of critical pressure in 
the corner. Therefore the second objective of these investigations is to analyse how the prediction of the 
developing corner flow changes, when a higher fidelity turbulence models that take anisotropy into account 
are applied. Therefore the second objective of these investigations is to analyse how the prediction of the 
corner flow changes for increasing incidence angles when the SA and the kcj-SST turbulence models with 
and without QCR extension as well as a Reynolds stress model are used. 



(a) Overview 


(b) Details, influence of grid resolution and topology, blue: 
negative c/, red: positive Cf 


Figure 4. Side-of-body flow separation at cruise design, Moo =0.85, 0 = 3.75°, Re = 5-10 6 . 



Figure 5. Influence of dissipation and turbulence models on CRM side-of-body flow separation, Mo© = 0.85, 
0 = 3.75°, Re = 5-10 6 . 


II. NASA Common Research Model and DPW-5 Test Cases 

For DPW-5 the NASA Common Research Model (CRM) civil transport aircraft configuration for cruise 
flight conditions (M^ = 0.85, Cl = 0.5, altitude H = 11300 m) is used as the reference geometry. The 
CRM optionally has a horizontal stabilizer as well as engines and pylons. In DPW-5 only the wing-body 
configuration is used as presented in figure [fl The CRM was designed by NASA’s Subsonic Fixed Wing 
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Technical Working Group and by Vassberg et al. 16 The wing has a slightly stronger pressure recovery at the 
last 10-15% local chord on the upper surface of the outboard wing section. The objective of this feature is 
to reduce boundary layer strength to control the development of a trailing edge separation and to create a 
challenge for turbulence models. The main geometrical features of the CRM are listed in table Tj Further 
details are published by Vassberg. 16 The geometrical and experimental data of the model can be found on 
the NASA CRM web site. 22 


Table 1. CRM geometrical data. 
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A. DPW-5 Cases 

DPW-5 required two mandatory test cases as a minimum. Additional parameter, grid, and turbulence model 
variations were allowed. 

1. Case 1, Common Grid Study: 

• Flow conditions: = 0.85, Cl = 0.5 ± 0.0001, Re = 5 • 10 6 

• Grids: sequence of common grids, custom grids 

2. Case 2, Buffet Study: 

• Flow conditions: M ^ = 0.85, Re = 5 • 10 6 , steady flow simulations 

• a = [2.5°, 2.75°, 3.0°, 3.25°, 3.5°, 3.75°, 4.0°] 

• Grids: medium (L3) common grid, custom grid 

III. Numerical Methods 


A. TAU CFD Solver 

Since the mid 1990s the Reynolds- averaged Navier-Stokes solver TAU is under development at DLR. It can 
be traced back to the German CFD project MEGAFLOW which integrated developments of DLR, aircraft 
industry, and universities. 23-25 Today the software package is under continuous development by the C 2 A 2 S 2 E 
department (Center for Computer Applications in AeroSpace Science and Engineering) of the institute and 
it is applied by DLR, European partners in industry and academia. 

TAU is an edge-based unstructured solver using the dual grid technique and fully exploits the advantages 
of hybrid grids. The numerical scheme is based on the Finite- Volume method and provides different spatial 
discretization schemes like central and upwind. 25 Here, a central scheme of second order accuracy, using the 
Jameson-type of artificial dissipation in scalar and matrix mode, has been applied. 26,27 Time integration has 
been performed using both, the explicit Runge-Kutta multistage and the Lower-Upper Symmetric Gauss- 
Seidel (LU-SGS) schemes. TAU has been developed with a particular focus on industrial aeronautical 
applications, thus providing techniques like overlapping grids for treating unsteady phenomena and complex 
geometries. Further details of TAU can be found in the reference. 25 
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B. Fluid-Structure-Coupled Simulation Procedure 

DLR’s fluid-structure-coupled (FSC) steady state simulation procedure, figure [S incorporates the in-house 
flow solver TAU and the computational structural mechanics (CSM) code NASTRAN® as main components. 
Additional modules included are a bi-directional interpolation routine for mapping aerodynamic loads to the 
structural nodes and transferring structural deflections back to the CFD mesh, and - closing the coupling 
loop - a volume mesh deformation algorithm. 



Figure 6. Simulation procedure for fluid-structure-coupled analyses. 

The analysis starts from an initial RANS CFD solution, which is computed on the undeformed grid. 
Then, static pressure and friction coefficients along with the identifiers, coordinates, and connectivity of the 
grid nodes, which constitute the CFD coupling surface, are transferred to the interpolation module using 
the Aerodynamic Mesh Interface Format (AMIF) specification. 

For each surface element in the CFD grid the interpolation module computes a force vector using pressure 
coefficient values, cell face area and cell orientation. Then, aerodynamic forces are mapped to the structural 
nodes located on the coupling surface. The corresponding finite-element surface data is provided from another 
AMIF file and processed in the same manner. Due to the considerable resolution difference, which usually 
exists between CFD and structural meshes, or when, as in this case, connectivity data of the finite-element 
surface nodes is not available, the application of a simple linear interpolation strategy is not applicable and 
a nearest neighbor search algorithm is used instead . 28 An assessment of both interpolation methods with 
respect to the coupling of aerodynamic forces between CFD and structural meshes is provided in . 29 For a 
given CFD face centroid i the nearest neighboring CSM grid point j is identified and a force component 
Fj,csM and associated moment Mj^csm = F%,cfd x are mapped to node j, figure ff\ This procedure 
ensures a conservative interpolation with respect to both force and moment balance on CFD and CSM side. 
An example showing a computed ^-distribution and the equivalent structural force distribution is given in 
figure El (a). 




Figure 7. Force mapping between CFD and CSM meshes. 
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(a) Interpolation of static surface pressure (top) to nodal 
forces (bottom) 


(b) Interpolation of structural deflections (bottom) to CFD 
surface mesh (top) 


Figure 8. Interpolation of aerodynamic loads and structural deflections. 


Next, nodal loads from the interpolation routine are re-formatted into NASTRAN® force cards and 
linked to the bulk data file. A linear, static structural analysis is performed and the resulting nodal deflection 
components along the coupling surface are mapped back to the CFD surface mesh, figure El (b). Because 
the nearest neighbor search algorithm used before is not appropriate for deformation fields, an interpolation 
scheme based on radial basis functions (RBF) is used . 28 The technique is particularly well suited for smooth 
functions , 30,31 like the deformations of aerodynamic structures considered in this application. 

Before a new flow solution is started, the interpolated surface nodal deflections are extrapolated into 
the volume mesh. This is achieved by applying the RBF interpolation functions used for the surface mesh 
deformation to the volume mesh nodes also. Additionally, the resulting deflections are superimposed with 
a weighting function based on wall distance in order to achieve a gradual decline of nodal deflections from 
the coupling surface into the flow field and to let them vanish for a specified distance, for example along the 
farfield boundaries. The method is applicable to both hybrid unstructured and block-structured meshes. 

Finally, a new CFD solution is computed on the deformed mesh. A typical convergence history for a fluid- 
structure-coupled simulation is plotted in figure gg The individual coupling steps are easily identified by steep 
increases in density residual and altered lift, drag, and pitching moment coefficient values. Iteration proceeds 
until user-defined convergence criteria, based on either flow or structural parameters, are accomplished. 



Iteration Step 

Figure 9. Coupled simulation convergence history. 
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DLR’s fluid-structure-coupled simulation approach has been validated using a variety of test cases and 
flow conditions, including both wind tunnel 32 and flight test data. 33 

C. Turbulence Modeling 

The DLR TAU code offers various turbulence models, ranging from simple eddy- viscosity to differential 
Reynolds stress models. The DLR results presented at the fifth AIAA Drag Prediction Workshop have been 
obtained using standard eddy viscosity models, i. e. the Spalart-Allmaras (SA) model 34 without the so-called 
fti and ft 2 term (SA-noft) and the ko;-SST model by Menter in its 1994 version. 35 

Such linear eddy viscosity models cannot predict the anisotropy of turbulent normal stresses near walls 
which is considered responsible for secondary corner flow phenomena. Indeed, using the SA and the kcj-SST 
model, a separation in the wing-fuselage intersection of the NASA CRM configuration is predicted at high 
incidence, which is not observed in the experiment. Improvement is expected by either non-linear exten- 
sions of eddy- viscosity models, e. g. explicit algebraic Reynolds stress models (EARSM), or full differential 
Reynolds stress models. In the following chapters the Quadratic Constitution Relation (QCR) extension in 
combination with the SA and the kcj-SST models as well as the DLR Reynolds stress turbulence model are 
described. 


1. QCR- extension of eddy -viscosity models 

As has been shown in DPW by Yamamoto et al. and Sclafani et ah, the corner flow separation can be 
suppressed, using the so-called Quadratic Constitution Relation (QCR) extension. 36-38 This modification 
has been presented by Spalart 39 as an example of how linear eddy viscosity models can be improved for 
predicting normal stress anisotropies. In general form it states the components of the Reynolds stress tensor 
as 

(i) 


75R _ -p( EVM) p( QCR ) 

pH Z j — pi i - • + P-^ij 5 


where pR\f VM ^ are the Reynolds stress components according to any linear eddy viscosity model and 
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is the QCR extension. The term 


Oij — 


dUj 

dXn 


dUj 

dXn 


/ 


dUrndUrr 
dx n dx n 


(3) 


denotes the components of the normalised rotation tensor with Vi representing the components of the mean 
velocity. The coefficient has been set to c n \\ — 0.3. 

Note that the above formulation depends on the respective eddy viscosity model, since the Boussinesq 
hypothesis reads 

pR% VM) = -2 + \pk5ij (4) 

where is the eddy viscosity, 


Sa 


c* 

Oij 



$ij ^ Skk&ij 


(5) 

(6) 


are the components of the simple and the traceless strain rate tensor and 5ij is the Kronecker symbol. The 
crucial point is the specific kinetic turbulence energy &, which is not available, e. g. with the SA model. For 
this reason the QCR extension {2) is cast into the following form 


Mf CR) = 4 C nllP (t) 


^ik^jk T Qjk^ik 


^mn^rrin T D mn D mn 


( 7 ) 
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which only contains the eddy viscosity ancHs independent of the availability of k. 

Furthermore, the production term of the ^-equation, employed e. g. by the ko;-SST model, deserves 
attention, as it reads 

pP (Z) = -pRkiS k i = - (pR[ f VM) + pR^ CR) ) S k i (8) 

As can be shown, 

pR$ CR) S kl = 0 (9) 

so that the QCR extension does not alter the transport equations of the underlying eddy viscosity model. 
Thus, the QCR extension can be easily implemented in a general form by adding the corresponding fluxes 
d/dxkipR^ 0 ^) and d / dxk(j)R^ CIV> Ui) to the momentum and total energy equation, respectively. 

The effect of the coefficient c n n on the predicted normal stress anisotropy can be analysed using standard 
assumptions on boundary layers. Assuming the wall-normal derivative of the wall parallel velocity component 
dU /dy being the only one non-zero component of the velocity gradient tensor, there are only two non-zero 
QCR stress components, i. e. 


?rT R) = 

dy 

pgg™ = 

dy 


( 10 ) 

( 11 ) 


where index 1 denotes the wall-parallel mean flow direction and index 2 the wall-normal direction. In the 
log- law 

( 12 ) 

(13) 


hypothesis, 

so that one finally obtains the QCR stresses 



= Kpu r y, 

dU 

u T 

dy 

Ky 

is the friction velocity. ' 

ku \ 

= = °- 09 


pR[® CR) = 2c nU fc, 
pR^ CR) = -2 c nll k. 


Combination with the Boussinesq hypothesis (3} yields 


pRn = 2k + Cniiy/c^j , 

PR22 = 2k f— + C n ii^/c^J 


(14) 

(15) 

(16) 

(17) 

(18) 


from which the corresponding components of the Reynolds stress anisotropy tensor 


are directly obtained as 



b i:i = Tl R - -Sij 
3 2k 3 

(19) 

bn 

= = 0.09, 

(20) 

b22 

CrZl \AV 0.09, 

(21) 

^33 

= 0. 

(22) 


As shown in Figure [101 computations for the flat plate, using the ko;-SST model with QCR extension, confirm 
the above analysis. Note that for the flat plate neither the skin friction nor the velocity profile is affected by 
the QCR extension. In contrast a small upstream shift of the shock is found for the flow around the RAE 
2822 airfoil, Case 9 and Case 10, 40 confirming findings by Yamamoto et al. 37 
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Figure 10. Flat plate, local Reynolds number R x = 9.15 • 10 6 . 
QCR extension of the kc^-SST model. 


Reynolds stress anisotropy components due to 


2. SSG/LRR-uo differential Reynolds stress model 

The SSG/LRR-o; model developed by DLR 41 is based on the Reynolds stress transport equation for com- 
pressible flow 


d 


{pRij ) 
dt 


+ (pRijUk^J — pPij + pHij — p€ij + pDij + pMij , 


where the production term is exactly given by 


- nP . __5 d El_R ^ 

P-Tij - n lk 


(23) 


(24) 


The Re-distribution term is modeled as 

/>n„ = -(c\pe+^C; Assa/LRR, mk}h J 


+C t 33 a/Lm^ ) + ( c^SG/LRR, _ c .,< SSO /LR*) ^ 


— 7. O* 

ij 


+c { ssg/lrr: )_ k ^ bikS . k + b . kSik _ l bmnSmnS .fj + d SSG/LRR) pk (bikWjk + b jk w ik ) , (25) 


where the dissipation rate e is computed from the specific dissipation rate u according to 


e = 


(26) 


with = 0.09 and the specific kinetic turbulence energy is related to the trace of the specific Reynolds 
stress tensor by 

k=^Ru. (27) 

The coefficient values vary from the LRR-values near the wall to the SSG-values at the outer edge of the 
boundary layer according to 


C (SSG/LRR ) = F ^Q(LRR) + ^ ^ Q(SSG) 

(SSG/LRR) TJ, x pi*,(SSG) 

~ f-L-^1 )Gi , 


(28) 


where F\ is Menter’s blending function. 35 The bounding values are listed in table 2. Note that the value of 
the LRR-parameter C ( LRR ^ = 0.52 has been recently modified. 42 


10 of ,26] 


American Institute of Aeronautics and Astronautics 


Table 2. Values of closure coefficients for the SSG and the LRR contributions to the SSG/LRR-id re-distribution 
term. cf RR) = 0.52. 



c\ 

C{ 

c 2 

d 3 

a 3 * 

a 4 

c 5 

SSG 

3.4 

1.8 

4.2 

0.8 

1.3 

1.25 

0.4 

LRR 

3.6 

0 

0 

0.8 

0 

18C (lrr) + 12 

-UC^ LRR) 


11 11 


The dissipation term is modeled as anisotropic tensor with components 


P e <;i ~ o e< +- 


Different from, 41 the diffusion term is modeled here by simple gradient diffusion 


pDij — _ 


d_ 

dx k 


jj | jj{SSG/LRR)Pj^ j &Rij 

+ / dx k 


(29) 


(30) 


for enhancing numerical robustness. In this, // is the averaged molecular viscosity, and the diffusion coefficient 
varies between the bounding LRR and SSG values according to 


2 C„ 


D (SSG/LRR) = Fia * + (1 _ Fi) ± ^ 

3 C/y 


(31) 


with cr* = 0.5 and C s m 0.22. 

The specific dissipation rate is provided by Menter’s baseline ca-equation, 35 where uj is limited similar to 
Durbin’s suggestion for eddy-viscosity models. 43 

According to Wilcox, 44 in the LRR part the following components of the anisotropy tensor, are obtained 

« 0.120 (32) 

« -0.114 (33) 

« -0.005, (34) 

where the little difference to zero trace of the numerical values is due to round-off errors. As one can see, 
the absolute values are larger than those associated with the QCR-extension, and in particular ^33 ^ 0. 

IV. Computational Grids 


&n 

^22 

^33 


8+12 C { 2 LRR) 

~ 33 c[ LR ^ 

_ 2 — 30 C^ LRR) 


A. CFD Grids 

1. Common Grids 

A six level common grids family of point-matched 0-0 topology multi-block grids has been build by Boeing. 45 
The sequence is based on an extra-fine grid (L5) with 40.9 • 10 6 hexahedral elements. To limit grid sizes a 
2-to-3 cell grid generation strategy was applied. The L6 grid (ultra fine) has been generated by refining L5 
by 1.5 in each parameter direction. The coarser grids L4-L2 and L3-L1 have been defined by dividing the 
grids by 8 appropriate for multigrid. 45 The derived medium grid (L3) with 5.1 • 10 6 elements represents a 
current grid size in industry for wing-fuselage configurations and will be used for the fluid-structure-coupled 
calculations and is presented in figure UR 
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2. DLR Custom Grids 



Since DPW-3 DLR is investigating un- 
structured prismatic and hexahedral el- 
ements dominated grids regarding their 
capabilities to be used for accurate drag 
predictions for aircraft 11, 12 with the soft- 
ware packages Centaur 46 and Solar. 47 
The hexahedral based approach offers po- 
tentially higher stretched elements, e.g. 
in wing span wise direction, because 
discretization errors are larger for pris- 
matic elements when highly stretched. In 
DPW-5 it was one objective of DLR to 
further compare the grid convergence be- 
haviour of the common grids with a hex- 
ahedral grid family generated with Solar. 

Because of the difficulty of adequate element shapes and sizes in corners with the Solar hex-dominated tech- 
nique, additionally three levels of Solar plus an overlapping full H-topology hexahedral grid block have been 
generated too. Due to the fact that Centaur is a very mature software, that prisms dominated grids are 
well established at the institute, and because Centaur offers hexahedral elements near the walls and partly 
in the field another objective was to compare the results of all four grid types (common hexahedral, Solar, 

TM 

Solar plus overlapping hex-block, Centaur ) with/without hexahedral wake blocks for Case 2. Results have 
been published in the frame of the the German Aerodynamics Workshop STAB. 13 Figure [K] shows the Solar, 

TM 

Solar plus overlapping, and Centaur grids. Here only the common grids L3, L4, and the Centaur grids 
with hexahedral wake block are applied. Grid sizes as well as the turbulence models used with these grids 
are listed in table [U 


Figure 11. Common hexahedral grid of DPW-5. 




(b) Solar + hexahedral overlapping (c) Centaur 

block 



(a) Solar 


Figure 12. Overview of DLR custom grids applied in DPW-5. 


Table 3. Grids and turbulence models applied for case 2, CFD grid nodes in million. 


Level 

Common 


Centaur 

+hexa-wake 


Turb 

Nodes 

Turb 

Nodes 

3 

4 

SA, SA+QCR, ku- SST, fcw-SST+QCR, SSG/LLR-w 
SA, SA+QCR, ku-SST, fcw-SST+QCR, SSG/LLR-w 

5.2 

17.4 

SA, SA-QCR 

37.4 


B. CSM Model 

A NASTRAN® finite-element structural model of wing, fuselage, horizontal tail plane, engine nacelles, and 
balance interface was made available by NASA Langley, figure [JJJ The model includes both right and 
left sides to account for the wind tunnel model’s non-symmetric inner structure. Joints between individual 
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components are modeled with rigid body elements. A rigid suspension is assumed at the balance interference. 
The finite-element discretization consists of 1.37 • 10 6 nodes, 6.82 • 10 6 elements, and 20.46 • 10 6 degrees-of- 
freedom. For the coupled simulations the engine nacelles and pylons were removed to more accurately 
represent the actual wind tunnel configuration. Coupling of aerodynamic loads between CFD simulation 
and finite-element analysis is established on the wing upper and lower surfaces. 



Figure 13. CRM finite-element model. 


V. Results 


A. FSC Simulations 

The purpose of the fluid-structure-coupled simulations here is to determine the static aeroelastic equilibrium 
state for one selected DPW-5 test case and to assess the influence of wing deformations on static pressure 
distributions and overall aerodynamic coefficients, in particular the pitching moment. This will help to 
quantify the effects, which have caused the observed deviations between CFD simulations and experiments. 

FSC simulations were run for DPW-5 test case 2 (cf. Chapter Aj) using the medium (L3) common grid and 
SA turbulence model. For an improved evaluation along with results obtained during recent investigations 
of support system effects, 20,21 the angle-of-attack is varied from 0.0 to 4.0°. 

Generally, the aeroelastic effects observed on the CRM were found to be larger than for the DLR-F6 
wing-body configuration 32 used during the Second and Third Drag Prediction Workshops. In figure 14’, 
the overall aerodynamic coefficients Cl and Cd obtained from the conventional CFD and FSC simulations, 
respectively, together with experimental data from the NTF wind tunnel test campaign (Test 197, Run 44), 
are plotted as a function of angle-of-attack. The coupled simulation results show lower overall lift coefficient 
values compared to the conventional CFD results, with the difference between FSC and CFD increasing 
with angle-of-attack. The lift reduction is due to the nose-down wing twist deformation induced by the 
geometric bending-torsion-coupling of the backward- swept wing as the wing is bent upwards by the external 
aerodynamic loads. 

The drop in lift starting between a = 3.0° and a = 3.5° found in the numerical data is due to an 
over-prediction of the side-of-body flow separation size in the medium (L3) hexahedral grid when using the 
SA turbulence model. The separation effect and related numerical issues were investigated and discussed 
previously in section El In the linear region, i.e. for a < 3.0°, deviations between the coupled simulation 
results and experimental data are considerably smaller than for the conventional CFD analysis, table® The 
remaining error is very similar to the differences found by Rivers et al. 20, 21 to be caused by the model support 
system. This suggests that including both aeroelastic and support system effects in the numerical simulation, 
together with a physically correct turbulence model, will allow for removing most of the previously observed 
deviations. 

Differences in drag coefficient between conventional CFD and FSC remain very small for incidence angles 
up to 2.5°. The deviation at off-design conditions is caused by the shock-induced flow separation on the 
outboard wing. In the coupled simulation, the separation starts to develop later, i.e. at a higher angle-of- 
attack, and, compared to the conventional CFD analysis, extends over a smaller spanwise portion of the 
outboard wing. This is due to the lower local angles-of-attack in the outer region of the deformed wing. 

For pitching moment coefficients, figure |15j the deviations between numerical and experimental data 
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Figure 14. Overall aerodynamic lift and drag coefficients for DPW-5 test case 2. 


Table 4. Lift coefficient deviations between numerical predictions and experimental results. 


a 

ACl,cfd 

ACl,fsc 

0 

0 

0 

0.0456 

0.0260 

0 

0 

cd 

0.0574 

0.0118 


is greatly reduced by taking into account wing deformation in the coupled simulation. Still, considerable 
differences remain, even around the design point. Again, including support system effects appears likely to 
move the numerical predictions closer to the experimental data. 



Figure 15. Pitching moment coefficient for DPW-5 test case 2. 


Figure 1T6 shows a comparison of chordwise static pressure distributions between CFD and FSC simula- 
tions and wind tunnel test data taken from the NTF campaign for four different spanwise wing sections at 
a = 3.0°. At the innermost section, figure HB1 (a), where wing deformations are very small, cf. figure 17 (b), 
both numerical methods are in good agreement with each other and the measured pressure distribution. 
Although twist deformation in this section is only about —0.0115°, the shock location predicted by the FSC 
simulation lies somewhat closer to the wind tunnel data. The mid- wing section, figure \M\ (b), already shows 
some effects of aero-elastic deformation between the leading edge and about 75% chord, with a decreased 
rooftop pressure level, reduced pressure along most of the wing lower side, and, again, a more precise shock 
location. At 77 = —0.727, figure CE2J (c), wing twist has increased to 5 = —1.09° and the differences between 
conventional CFD and FSC become even more apparent. Here, only the coupled simulation is in good agree- 
ment with both measured rooftop pressure levels and shock location. At the outermost section, figure 16 (d), 
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the shock location as predicted by the FSC simulation has moved significantly downstream and a double- 
shock pattern has developed. Twist deformation has increased to e = —1.41° considerably reducing the 
local incidence angle. As a result, the pressure distribution at 77 = —0.950 resembles those for conventional 
CFD computations at lower angles-of-attack, cf. figures \2\ (b) and||(b), where a similar double-shock system 
exists. Unfortunately, the true shock position can not be determined from the experimental data due to an 
insufficient spatial resolution of pressure taps. 




(a) 77 = 0.131 


(b) 77 = 0.502 




(c) 77 = 0.727 


(d) 77 = 0.950 


Figure 16. Chordwise static pressure distribution at a = 3.0° for different spanwise wing sections; TAU only 
(CFD), fluid-structure-coupled (FSC), and NASA NTF experiments (NTF). 


In figure !J7] the wing bending and twist deformations are plotted as a function of angle-of- attack at wing 
tip (a) and as spanwise distribution for a = 3.0° (b). As previously seen with lift coefficient, figure |14 (a), 
good linearity exists for a < 3.0°. Between a = 3.0° and a = 3.5° the onset of the side-of-body flow 
separation can be identified by small decreases in both bending and twist deformations. The declining slope 
for a > 3.0° is caused by the growing shock-induced flow separation on the outer wing. 

The spanwise wing twist and bending distribution for a = 3.0° is plotted in figure [17] (b). Due to the 
fact that from 77 ~ 0.40 outward the c/4- line lies behind the model reference center, any aeroelastic wing 
deformation will not only result in a change of spanwise lift distribution, but also reduce the overall nose- 
down pitching moment. Unfortunately, no experimental deformation data was available for comparison as 
the corresponding wind tunnel test was still ongoing at the time of publication of this paper. 


B. Turbulence model study 

In this section the influence of turbulence models on the prediction of side-of-body flow separation for the 
CRM wing-body configuration is investigated. Here two eddy viscosity models, i.e the one-equation model 
by Spalart-Allmaras (SA) and the shear stress transport (ko;-SST) k — u two-equation model by Menter, 
and the differential Reynolds stress model (SSG/LRR-cj) developed at DLR are applied. For this study the 
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(b) Spanwise wing deformation 


Figure 17. Wing Bending and Twist Deformations. 


CRM hybrid grids of medium level (L3) with 5 million nodes and of fine level (L4) with 17 million nodes 
have been used. On the medium grid, computations are performed with the SA, SA with QCR extension 
(SA+QCR); kcj-SST, kcj-SST with QCR extension (SST+QCR); and the SSG/LRR-o; models, whereas the 
predictions on the fine grid are obtained by the SA, SA+QCR and the SSG/LRR-cj models only. In addition 
the SA and SA+QCR models are also applied in combination with the Centaur grid to verify the effects on 
a prismatic element dominated grid too. 

In this study, the viscous fluxes of main and turbulent equations are discretized using central differences. 
The inviscid fluxes of the main equations are calculated by a central scheme with matrix dissipation. For 
the turbulent equations the convective fluxes are approximated with a central scheme for the SA model and 
with second-order Roe’s scheme for the kcj-SST and SSG/LRR-cj turbulence models. Steady computations 
are performed using a semi - implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method. For the 
computations here, slightly different numerical dissipation parameters than typical TAU dissipation settings 
have been used. The parameters used here lead to slightly more dissipation. This is done in-order to be able 
to obtain a steady converged solution with all three turbulence models on the grids employed. 

Simulations were performed for the CRM wing-body configurations from a = 2.5° to a = 4.0° and for 
selected cases up to a = 4.25° with incidence angle increment of Aa = 0.25°. For the cruise design condition 
(CL=0.5) simulation is carried out using the solution of a = 2.5°. 

In figure 18, the force and moment coefficients predicted by different turbulence models on the L3 grid 
are shown. In the figure experimental data are used for reference purposes. With all the turbulence models 
continuous increase in lift as the incidence angle increased is predicted - no lift breakdown is observed. The 
SA model predicted higher lift at all the incidence angles than the kcj-SST and SSG/LRR-o; models. Below 
the incidence angle of 3.5° the kc^-SST model predicted lower lift than SSG/LRR-o;. For a > 3.5° the 
ko;-SST and SSG/LRR-o; delivered almost the same lift coefficient. The eddy viscosity models with QCR 
extension delivered lower lift compared to the corresponding non- QCR models. 

The drag force polar demonstrates that the SA model delivers lower drag than the kcj-SST and SSG/LRR- 
cj models and the drag predicted by the latter model is in between the predictions delivered by the kcj-SST 
and SA models. With QCR extension both the SA and kcj-SST models predict higher drag than the 
corresponding models without QCR extension. 

The pitching moment coefficient is displayed in figure 18(c) for different turbulence models. The SA 
model predicts the lowest pitching moment coefficient compared to the kc^-SST and SSG/LRR-o; models. 
The predictions delivered by SSG/LLR-cj can be found in between the predictions of the SA and kcj-SST. 
The eddy viscosity models with QCR extension predicted a higher pitching moment coefficient compared to 
non- QCR models. 

In figure [19 the Cp distribution predicted on the L3 grid at a = 3° and a = 4° is compared for different 
turbulence models at two spanwise sections of the wing. Here the experimental data are used as reference. 
Differences persist in the location of the shock and in the pressure distribution in the separated flow region. At 
a = 3° the SA model predicts the location of the shock slightly downstream of the SSG/LRR-cj predictions. 
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Figure 18. Comparison of force and pitching moment coefficients for different turbulence models on L3 and 
Centaur grid. 


The kcj-SST model predictions unveil that the shock is predicted slightly upstream of the shock location 
delivered by the SSG/LRR-cj. The upstream shift of the shock by the kcj-SST model is about half of the 
difference between the SA and SSG/LRR-Cc; models. At a = 4° the SA model predicts lower pressure on 
the suction side of wing than the other turbulence models and delivers the shock slightly downstream of the 
kcj-SST and SSG/LRR-cj predictions. At this incidence angle very small difference between the kc^-SST and 
SSG/LRR-cj with regard to shock location is observed. 

With QCR extension the SA and ka;-SST models shift the shock location slightly upstream compared 
to the corresponding non-QCR models. This upstream movement of the shock caused larger shock induced 
separation on the main and outboard of the wing which lead to lower lift and higher pitching moment (see 
figure T8Tc)). 

To demonstrate the influence of grid refinement on the predictions, drag and moment coefficients predicted 
on L3 and L4 grids are shown in figure [201 for the SA, SA+QCR and SSG/LRR-cj models. The predictions 
obtained by the SSG/LRR-o; unveil grid independence up to the incidence angle of 3.5°. On the L4 grid the 
SA model with and without QCR predicted slightly higher lift and lower pitching moment coefficients than 
on L3 grid. 

In figure HI] Cp distributions predicted on the two different grids by the SA+QCR and SSG/LRR-cj models 
at a = 3° are shown here. The SA model predictions are not depicted as the difference between predictions 
on both grids is of same magnitude that was observed with SA+QCR and hence only the SA+QCR results 
are discussed. For the SSG/LRR-o; Reynolds stress model no significant influence of grid refinement was 
observed, whereas with the other model downstream movement of shock was predicted on the fine mesh. 
This is the possible reason for slightly higher lift and lower pitching moment coefficients on the L4 mesh by 
the SA and SA+QCR models (see figure [20]) . Please note that no fluid-structure coupling is applied for the 
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x/c 

(c) Tj = 0.50 (a = 4°) 


x/c 

(d) 77 = 0.95 ( a = 4°) 


Figure 19. Cp distribution at spanwise sections of the wing for different turbulence models on L3 grid at a = 3° 
and a = 4°. 


turbulence models investigation. As observed in DPW a large discrepancy of experimental and numerical 
data exists at 77 = 0.95 due to reasons investigated by Rivers et al . 20,21 

To examine the influence of turbulence models on the side-of-body flow separation surface skin friction 
lines for the a = 3° case are displayed in figure. \22\ On the L3 grid a very small separation is predicted with 
the SA model whereas with the SA+QCR, ko;-SST, kcj-SST+QCR and the SSG/LRR-cj models no separation 
in the side-of-body region is predicted. On the Centaur grid a sightly larger side.of.body separation can be 
observed for the SA model due to the significant higher grid resolution in that area compared to the L2 grid. 
Applying the SA+QCR model again a reduction of the size of the separation is achieved. As the incidence 
angle increases, the SA model predictions show an increased size of side-of-body flow separation. However, 
the size of flow separation is not as big as it is predicted by Yamamoto et al . 37 and Sclafani et al . 38 on 
the same hybrid L3 grid. The reason for this small separation is marginally higher dissipation which in the 
current investigation arose from slightly different numerical dissipation parameters that were different from 
typical TAU settings but that had been necessary to be consistent for all grids. Despite small separation 
in the side-of-body region comparisons are made for different turbulence models. Surface skin friction lines 
obtained at a = 4° by different turbulence models are shown in figure 231 The SA model predicted large 
separation whereas the kcj-SST and SSG/LRR-cj predicted no separation in the region. With QCR extension 
the side-of-body flow separation is reduced by the SA model. The kcj-SST+QCR model predicted a similar 
flow topology as the kcj-SST predictions. 

Predictions obtained on the L4 grid by the SA, SA+QCR and SSG/LRR-o; models are depicted in 
figure [Ml The SA model predicted larger flow separation on the L4 grid than on the L3 grid at a = 3°. 
At higher incidence angle much larger separation is predicted by the SA model. This separation extent 
is reduced by the QCR extension (see figures [25). The SSG/LRR-o; has not predicted any side-of-body 
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(d) Pitching moment (SA) 


Figure 20. Comparison of force and pitching moment coefficients for the SA and SSG/LRR-w models on L3 
and L4 grids. 


separation on the grid. To illustrate the size of the separation in the wing-body junction in comparison with 
the wing size a detailed view of the wing with skin friction lines is displayed in figure 26. 

In the experiments of the CRM no separation in the side-of-body region was observed. Such trend is re- 
produced with the SSG/LRR-cj. With QCR extension, which accounts for the anisotropy of normal Reynolds 
stresses, the separation extent is reduced but it is not completely eliminated in the present investigation. 

VI. Conclusions 

the influence of wing deformations on static pressure distributions and overall aerodynamic coefficients of 
NASA’s Common Research Model were investigated using the medium (L3) common grid, the SA turbulence 
model, and a finite-element structural model kindly provided by NASA Langley. Static fluid-structure- 
coupled simulations were run at M ^ = 0.85 and Re = 5-10 6 with the angle-of-attack varying between 0.0 
and 4.0°. Numerical results were compared to experimental data from a wind tunnel test campaign in 
NASA’s National Transonic Facility. Generally, the deviations of lift, drag, and pitching moment coefficients 
observed between DPW-4 and DPW-5 computational results and measured data are considerably reduced 
by taking into account elastic wing deformations. Lift coefficient values predicted by the coupled simulation 
are lower than for the conventional CFD computations, leading to considerably smaller deviations from 
the experimental data. For drag coefficients, significant differences between the conventional CFD and 
FSC analyses only occur at off-design flow conditions and are mostly due to variations in the development 
of the shock-induced separation on the outboard wing. Deviations between numerical and experimental 
pitching moment coefficients are substantially reduced by taking into account wing deformation. Due to 
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(b) r) = 0.95 (a = 3°) 
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Figure 21. Cp distribution at spanwise sections of the wing for different turbulence models on different grids 

at o = 3° and a = 4°. 


the pitching moment’s strong sensitivity with respect to the overall static pressure distribution, differences 
remain relatively large. Regarding chordwise static pressure distributions, some minor aeroelastic effects 
become visible in the mid-wing section, increasing in magnitude towards the wing tip. In general, the FSC 
simulations provide a significantly more accurate prediction of rooftop pressure levels, pressure distribution 
on the wing lower side, and shock location. Wing bending and twist deformations show a good linearity for 
a < 3.0°. For higher angles-of-attack, wing deformations are influenced by the side-of-body flow separation 
and the shock-induced flow separation on the outer wing. Based on the results found in this study and by 
Rivers et ah, 20,21 it is suggested that further numerical investigations should include both aeroelastic and 
support system effects, together with a high-quality turbulence model. 

In a second part, a turbulence model study was performed for the CRM wing-body configuration on 
medium (L3) and fine (L4) hybrid grids provided by DPW-5 committee. Here eddy viscosity models, i.e 
the SA and kcj-SST models, with and without QCR extension and the SSG/LRR-o; differential Reynolds 
stress model were applied. On the L3 grid the kcj-SST model predicted lower lift and higher pitching moment 
coefficient than the SA and SSG/LRR-cj models. The aerodynamic coefficients delivered by the SSG/LRR-cj 
were in between the prediction of the SA and kcj-SST models. With QCR extension, the SA and kcj-SST 
predicted the shock slightly upstream compared to the corresponding non- QCR models and hence larger 
shock induced separation on the main and outboard of the wing, which lead to lower lift and higher pitching 
momentum coefficients. The flow separation in the side-of-body region predicted by the SA model was very 
small at cruise design conditions which increased in size as the incidence angle increased. On the L4 grid 
much larger separation was predicted than on the L3 grid at the corresponding incidence angle. The QCR 
extension in combination with the SA model only reduced size of separation but it did not eradicate side-of- 
body flow separation. The force, pitching moment coefficients and the surface pressure distribution predicted 
by the SA model with and without QCR extension were found to be grid depended. On the L4 grid slightly 
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(a) SA (Centaur) 


(b) SA (L3) 


(c) SST (L3) 




Figure 22. Comparison of side-of-body separation at a = 3° for different turbulence models on L3 and Centaur 
grids. 
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Figure 23. Comparison of side-of-body separation at a = 4° for different turbulence models on L3 grid. 
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Figure 24. Comparison of side-of-body separation at a = 3° for different turbulence models on L4 grid. 





(a) SA (b) SA+QCR (c) RSM 


Figure 25. Comparison of side-of-body separation at a = 4° for different turbulence models on L4 grid. 
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Figure 26. Comparison of skin friction lines on the wing at a = 4.0° for different turbulence models on L4 grid. 
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higher lift, lower drag and moment coefficients were predicted and the shock was located slightly downstream 
compared to the predictions on L3 grid. With the kcj-SST model both with and without QCR extension 
side-of-body flow separation was not observed on the L3 grid. The lift, drag, pitching moment coefficients 
and the surface pressure distribution predicted by the SSG/LRR-cj model on both grids were observed to 
be grid independent. On both grids the SSG/LRR-cj predicted attached flow in side-of-body region which 
is in agreement with the experimental observations. The results on the DLR custom Centaur grid showed a 
slightly larger side-of-body separation can be computed with the SA model due to the higher grid resolution 
compared to the L3 and L4 grids. The application of the SA+QCR model reduces the separation, similar 
as shown on the common grids. 
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